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Abstract. We present the results of the BeppoSAX observations of two BL Lacs, OQ530 and S5 07 16-1-7 14, as part of a ToO 
program for the simultaneous observation at radio, optical. X-ray and TeV energies. Both sources are detected in the LEGS and 
MEGS, with S5 0716-1-714 visible also in the PDS band, up to about ~ 60 keV. The X-ray spectra of both sources are better fitted 
' by a double power-law model, with a steep soft X-ray component flattening at harder energies, with breaks at 0.3 and 1.5 keV, 

^■f^ ' respectively. The concave shape of the spectra in both objects is consistent with soft X-rays being produced by the synchrotron 

, and harder X-rays by the inverse Gompton processes. Also the X-ray variability properties confirm this scenario, in particular 

for S5 0716-1-714 our observation shows variations by about a factor 3 over one hour below ~ 3 keV and no variability above. 
Their simultaneous broad band energy spectral distributions can be successfully interpreted within the frame of a homogeneous 
synchrotron and inverse Gompton model, including a possible contribution from an external source of seed photons with the 



O ' different spectral states of S5 0716-1-714 being reproduced by changing the injected power. The resulting parameters are fully 



■ 

X 



consistent with the two sources being intermediate objects within the "sequence" scenario proposed for blazars. 
Key words. BL Lacetae objects: general - X-rays: galaxies - BL Lacetae objects: individual: OQ530, S5 0716-1-714 



^ 1 . Introduction spectral features: with respect to the presence or the absence 

d .... of broad emission lines (EW> 5 A) they are classified as Flat 

■ - Blazars are radio-loud Active Galactic Nuclei whose spectra c » n j- r\ m t u- . i r>T 

^ Spectrum Radio Quasars or as BL Lac objects respectively. BL 

are characterized by variable nonthermal continua emitted by ^ c i c j ittit t dt j- . 

■' ■' Lacs are further classmed as HBL or LBL according to the po- 
relativistic lets oriented close to the line of sight. They dis- . . , ^ i t nt i • ^- ^ T^\r^. j 
■' a J sition of the synchrotron peak: LBLs peak in optical-UV band 
play double peaked Spectral Energy Distributions (SED) ex- ... . „ , ■ o r^- ■ innc\ t,, 
^ f f bj \ ' while HBLs in X-rays (Padovani & Giommi, 1995). The ex- 
tending from radio up to 7-rays and sometimes to TeV fre- . 1 j c^, or-i^ j u- u • u v. 

° r / ^ tremely wide energy range of the SED and the high variability 

quencies. The first peak (radio to UV/X-rays) is generated by 11. i.- 1 u 

^ r V J J b ■'of these sources makes simultaneous multi wavelength observa- 

synchrotron emission produced by relativistic electrons flow- ^ w w j . j • u • rr, 

-' f J tions a powerful tool to understand their physics. Therefore we 

ing along the let, while the second component (extending to . , ^t, ^ c r\ . -.n cavu 

° ° f \ b proposed a number of Target of Opportunity BeppoSAX ob- 

7-rays) is commonly attributed to inverse Compton emission. f,, • u- u . . f v j- . j -.u 

' ■' ' ■> ^ servations of blazars in high state of activity, coordinated with 

The same synchrotron emitting electrons upscatter "seed" pho- 1 j • ■ u ■ • f.u 

■' b f f optical and TeV monitoring campaigns: the brightening of the 

tons to higher energies. The origin of this seed radiation field is . > j u.- n cavmt-t 

° ° ° source in one of these bands would trigger BeppoSAX NFI ac- 

still debated: it could be the synchrotron radiation itself (SSC . .^.u- ij .jn cav 

^ tion. As part of this program we already presented BeppoSAX 
models, Maraschi,Ghiselhni&Celotti, 1992) or produced out- , . ^ -m ni^c onnc /ion /t- if • . 1 

., , . , , , , ... observations of ON 231, PKS 2005-489 (TaghafeiTi et al. 

side the lets by different mechanisms (ERC, external radiation ^rvrvrv --irvrvix j nr t /i-. • .1 '^nn'^\ t .u- 

^ , „ „ l^^. ^, • „. • „ 2000, 2001) and BL Lac (Ravasio et al. 2002). In this paper 

Compton, see Sikora, Begelman & Rees, 1994; Ghiselhm & 1. f.u n cav u f<~.<~.con 

^ ' ' o we report the results of the BennoSAX observations of OQ 530 

Madau, 1996; Dermer & Schlickeiser, 1993; Blazeiowski et al., jccmi/c ti/i/.u c »u j cf^u » • c \ 

and S5 07 loH-7 14 (the fourth and hfth triggers of our program). 

2000). Blazars are divided in two subclasses according to their 

® OQ530 — It was optically identified in 1977 (Kuhr 
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NRAO-Bonn radio survey at 5 GHz. In 1978, it was classi- 
fied as a BL Lac by Miller (1978), because of its stellar aspect 
and the faintness of its optical spectral features. The redshift 
z=0.152 has been determined from absorption features of the 
host galaxy as well as from [Oil], A = 3727 A and [OIII], 
AA = 4959,5007A emission lines (Stickel et al. 1991;Stickel 
et al. 1993), confirming previous estimates by Hagen-Thorn 
& Marchenko (1989). The host galaxy of OQ530 was stud- 
ied by Abraham et al. (1991), Stickel et al. (1993), Wurtz et 
al. (1996) and Scarpa et al. (2000). The data given by these 
authors, transformed into Rc (Cousins) give sUghtly different 
values for the host galaxy magnitude: 16.8, 16.2, 16.6 and 16.1, 
respectively . The first three used ground based instruments and 
found the galaxy to be probably an SO, the last ones used HST 
observations and found the galaxy to be an elliptical, like all 
other known host galaxies of BL Lac sources. In our monitor- 
ing, started in 1994, we never observed OQ530 fainter than 
Rc = 15.8, so all these values are formally acceptable for the 
luminosity of the host. Given the systematic uncertainties of 
these estimates we adopt here an average value of Rc = 16.5. 
The effective radius of the galaxy is found by all authors to be 
less than 4", so that it is nearly fully included in our photomet- 
ric aperture (radius 5" ). Our BeppoSAX pointing was made 
when OQ530 was at Rc = 14.8, so the host contribution to 
the optical flux is not negligible: we estimated this contribution 
assuming the spectral shape and k-correction for an elliptical 
galaxy by Fukugita et al. (1995), i.e. B =18.59, V =17.17, 
Rc =16.50 and /c =15.74. 

In the X-ray band OQ530 has been detected by the 
Einstein Observatory in May 1980. The data are well fitted by 
a power law spectrum with a — 0.48lJ;3 in the 0.2-3.5 keV 
range and a flux Fikev = 0.2lAtJy (Worral & Wilkes 1990). 
Some years later, in May 1984, a three days EXOSAT obser- 
vation highlighted a non variable behaviour of the source on 
short time scales, with a flux of FikeV — 0.15 ± 0.02/iJy 
(Giommi et al. 1990). ROSAT observed this BL Lac object in 
July 1990. The 0.1-2.4 keV spectrum is well fitted by a sin- 
gle power law with an energy index a = 1.04 ± 0.05 and flux 
FikeV = 0.32fiiy (Comastri et al. 1995). BeppoSAX already 
observed OQ530 in February 1999. The source was detected 
by the LEGS and MEGS (< 10 keV) but not by the PDS. 
A single power law model fit the data (a — 0.55^032), but 
the residuals towards low energies suggest a concave spectrum 
(Giommi et al. 2002). 

The optical light curve of OQ530 from 1905 until 1977 
was derived by Miller (1978) using the Harvard College plate 
collection. A further study from 1966 to 1980 was performed 
by Barbieri et al. (1982) using the Asiago Schmidt plates. A 
reanalysis of these works, and an updating of the light-curve 
until 1997 was made by Nesci et al. (1997). A collection of 
published data until 1993 has been recently published by Fan & 
Lin (2000). The source shows an average monotonic decreas- 
ing trend since the beginning of 1900, when it was at B ~12, 
until 1995 when its average value was B ~16, with fast and 
irregular variations of about 1.5-2 magnitudes. It is still moni- 
tored by some of us (Perugia and Roma groups). On the night of 
February 14, 2000, we found the source very bright in the op- 
tical with Rc =14.3. Thus, we triggered the BeppoSAX TOO, 



but due to constraints on the satellite pointing, the source was 
observed only on March 3, 2000, when it had a magnitude of 
Rc =14.9. Moreover, due to problems with the satellite, the 
observation lasted less than 27 ks, instead of the 50 ks allo- 
cated. A second pointing was then performed on March 26, for 
about another 23 ks; during this second pointing the Rc magni- 
tude was 14.8. In both cases the source was weak in the X-rays 
and detected only up to 10 keV (see section 2). 

S5 0716+714 — It was discovered in 1979 as the opti- 
cal counterpart of an extragalactic radio source during another 
Bonn-NRAO survey (Kiihr et al. 1981). Two years later it was 
classified as a BL Lac by Biermann et al. (1981) because of its 
featureless optical spectrum and high linear polarisation. It is 
one of the brightest and most variable BL Lac objects, but, de- 
spite the many attempts to measure its redshift (the last one by 
Rector & Stoke (2001) with the Keck telescope), only a lower 
Umit z > 0.3 has been estimated by Schalinski et al. (1992) 
and Wagner et al. (1996) on the basis of the non detection of a 
host galaxy in deep images. 

S5 0716+714 was first detected in the X-ray band by the 
Einstein Observatory, but the flux was too low for a detailed 
spectral analysis (Biermann et al. 1981). In 1991 March 8- 
11, the source was observed several times with the ROSAT 
PSPG detector (0.1-2.4 keV) (Gappi et al. 1994; Urry et al. 
1996; Wagner et al. 1996) and was found to be bright and 
rapidly variable with a mean flux of ~ 10~^^ erg cm~^ s~^. 
S5 0716+714 has been already observed twice by BeppoSAX, 
in 1996 November 14 and in 1998 November 7 (Giommi et al. 
1999) but the source was not very bright and it was detected 
only up to 10 keV. The BeppoSAX SEDs showed the presence 
of a steep power law component below ~ 2.5 keV with a spec- 
tral index ai ~ 1.5 becoming harder towards higher frequen- 
cies, with an index a2 ^ 0.85. The former component was also 
found to be variable in correlation with the optical flux on a 
time scale of a few hours (Gionmii et al. 1999). 

S5 0716+714 has been regularly monitored in the optical 
since the end of 1994; a very bright level was recorded in 
February 1995 (Ghisellini et al. 1997), when it reached the 
magnitude of 12.78 in the Rc bandpass, while in the fol- 
lowing two years it has never been observed brighter than 
Rc = 13.0. A new strong outburst occurred in September 
1997 when S5 0716+714 was even brighter than in 1995 at 
i?c=12.58 (Massaro et al. 1999). After a relatively quiet pe- 
riod, in October 2000 we measured an optical flux compara- 
ble to that of September 1997 and a BeppoSAX ToO was ac- 
tivated. Unfortunately, for technical reasons the pointing was 
performed a week after the occurrence of the optical maximum 
and the brightness of S5 0716+714 had already decreased by 
about half a magnitude. In any case, it was brighter than in 
the two previous BeppoSAX observations, when it was around 
i?c=13.8, and, for the first time, it was detected in the hard 
X-rays by flie PDS up to about 60 keV. 

The paper is organised as follows: in section 2 we present 
the BeppoSAX observations and data analysis and in section 3 
the radio, optical and TeV observations; in section 4 we present 
the SEDs of the two sources in the framework of an SSC model. 
The discussion is given in section 5. 
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2. BeppoSAX Observations and Data Reduction 

BeppoSAX observed OQ530 (1418+546) twice during 2000 
(March 3^, March 26-27). In both observations the source 
was not detected by the PDS. 

S5 07 16+5 14 was observed for about one day from October 
30 to October 31, 2000; and it was detected up to 60 keV. 

Standard procedures and selection criteria were applied to 
the data to avoid the South Atlantic Anomaly, solar, bright 
Earth and particle contamination using the SAXDAS v. 2.0.0 
package. Exposures and count rates of the observations are re- 
ported in Table [2 

Data analysis was performed using the software available 
in the HEAsoft Package (XIMAGE 3.01, XSELECT v2.0, 
XSPEC 11.0.1, XRONOS 5.16). LEGS and MECS images dis- 
played a bright pointlike source: events for spectral and timing 
analysis were extracted from circular regions centered on the 
source, with radii of 8 and 4 arcmin, respectively. We extracted 
LEGS and MEGS background events in the same way, from 
off source regions of the field of view and checked the con- 
stance of the count rates during the whole observation. In any 
case, since LEGS and MEGS backgrounds are not uniformly 
distributed across the detectors, we choose to use background 
files extracted from long blank field exposures, available from 
the SDG public ftp site (see Fiore et al., 1999). The spectral 
analysis were performed with the XSPEG 11.0.1 package, us- 
ing the updated (01/2000) response matrices. 

The PDS was operated in the customary collimator rock- 
ing mode, where half the collimator points at the source and 
half at the background and they are switched every 96 s. The 
background-subtracted PDS spectrum for S5 0716+714 was 
obtained from the standard pipeline analysis carried out at the 
BeppoSAX Science Data Genter. 

2.1. OQ530 

Because of the faintness of the source and the short duration 
of the runs the source was not detected by the PDS, as can be 
seen from Table Q and the Scr upper limit does not add any 
useful information to the X-ray spectrum (e.g. for the SED in 
Fig. 6 below). From the comparison of LEGS and MEGS data 
with background files we decided to constrain further the en- 
ergy range of our analysis to 0.13-8.7 and 0.13-10 keV for the 
3^ and 26-27 of March observations, respectively. 

For both observations we have poor statistics at low ener- 
gies (LEGS band) and it is difficult to obtain a reliable fit to the 
observed spectrum. We rebinned the LEGS energy channels in 
order to have more than 25 net counts in each bin: this is a 
good compromise between having a sufficient number of bins 
and significant number of counts in each of them. The MEGS 
detected a slightly higher number of events, so we rebinned 
them to have 30 net counts in each bin. 

We fitted our spectra adopting the galactic interstellar ab- 
sorption TVh = 1.18 X 10^" cm^^, as measured by Dickey & 
Lockman (1990). Fitting the 3-4 March spectrum with a sin- 
gle power law model leaves large residuals at low energies and 
the quality of this fit is evidently poor (Xr/d.o.f. = 1.64/23). 
Therefore we tried to fit the spectrum with the sum of two 



power law models and we obtained a very steep soft X-ray 
spectral index ai — 6.7l" g and a flatter one towards higher 
energies a2 — 0.40 ± 0.2 (throughout the paper the errors are 
at 90% confidence interval for one parameter of interest, i.e. 
Ax^ = 2.7). The two power laws cross at Eb ^ 0.3 keV. This 
model gives a better fit {y^/d.o.f. = 0.96/21, F-test probabil- 
ity > 99.9%), but because the break is near to the soft limit of 
our energy range and we have only a few counts in this part of 
the spectrum, we can not really constrain the soft component of 
the two power law model. This soft X-ray excess is probably 
due to the very steep tail of the synchrotron component. At en- 
ergy higher than ~ 0.5 kcV, instead, we are probably observing 
an inverse Gompton component. 

We repeated the procedure for the observation of 26-27 
March. A single power law fits the data better than in the 
first observation, however towards low energies hints of pos- 
itive residuals are still evident. We obtain a — 0.75 ± 0.20 
(Xr/d-O-f- — 1.03/21). The measured fluxes between the two 
observations are quite similar, -F2-iofcey ~ 10~^^ erg cm~^ 
s^^ (see Table |2j. We performed the fit also with the sum of 
two power laws. As in the previous observation, in the soft X- 
ray band we find a very steep component, with practically the 
same spectral index ai — 6.7, but with even larger uncertain- 
ties (basically we can not constrain this component). At higher 
energies the second spectral index is a2 = 0.65 ±0.2. With this 
model we obtain a lower Xr (Xr = 0.95), but an F-test reveals 
that the improvement of the fit in this case is only marginally 
significant (^ 80%). In Table |2] we report a summary of our 
analysis. 

Also during the 1999 BeppoSAiX observation, Giommi et 
al. (2002) have an indication of a better fit with a broken power 
law (ai — 1.3, a2 = 0.4, and Ei, = 1.8 keV). However, their 
statistics is even lower and the fit improvement is marginally 
significant. In Table |3l we report all the published spectral 
parameters for the OQ530 X-ray observations. The spectral 
shape of the source seems to be almost constant, except for 
the ROSAT observation, that found a steeper spectrum. This is 
probably due to a mixture of the soft and hard components, that 
the PSPG was not able to disentangle. The 1 keV flux density is 
quite variable; the faintest measured 1 keV flux (by BeppoSAX, 
in 1999) is in fact more than three times weaker than the highest 
one (ROSAT 1990). 

We performed the temporal analysis of our observations 
in the two bands 0.2-1.8 keV LEGS and 1.8-6 keV MEGS. 
Because of the faintness of the source we choose a time bin- 
ning of 1 hour, in order to reduce uncertainties caused by the 
poor statistics. During both observations we did not detect vari- 
ability. 

2.2. S5 071 6+714 

During the BeppoSAX observation of October 2000, 
S5 0716+514 was detected also by the PDS, so we could 
perform our analysis on a wider energy range (0.15-57 keV). 

We performed a preliminary analysis on LEGS+MEGS and 
MEGS+PDS data alone, in order to get some useful hints 
for fitting the whole spectrum. Initially we fitted the spec- 
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OQ530 





LECS 


MECS 


PDS 


Date 


exposure (s) count rate'' 


exposure (s) count rate*^ 


exposure (s) count rate'^ 


3-4 March 2000 


19625 (7.4 ± 0.8) x 10"^ 


26633 (1.3 ±0.1) X lO-"" 


10270 (1 ± 6) X 10"^ 


26-27 March 2000 


17181 (8.2 ± 0.9) X lO""* 


22718 (1.1 ± 0.1) X 10"^ 


10361 (2±6)xl0"^ 



S5 0716+714 



30-31 October 2000 


19192 (4.7±0.2) X 10~^ 


43459 (5.6 ± 0.1) x 10~^ 


22444 (0.12 ±0.04) 



Table 1. Journal of BeppoSAX observations. ^ 0.1-10 keV; ^ 1.5-10 keV; 12-100 keV. For OQ530 the PDS count rates are 
lower than the uncertainties and we excluded these data from our analysis. 



OQ530 



Date 


Ql 




02 


(xlO"'') 


FikeV MJy 
(1) (2) 


TP** 

(xlO-^2) 


Xi/d.o.f. 


3-4 March 


0.55 ±0.2 


9 1+U.fi 
^•-■--0.5 






0.14 


1.1 


1.64/23 


3^ March 


' -1.8 


~ 10~^ 


0.40 ± 0.2 


1.8 ±0.5 


~ 10~* 0.12 


1.2 


0.96/21 


26-27 March 


0.75 ±0.20 


9 C+U.5 






0.16 


0.9 


1.03/21 


26-27 March 


6.7±? 


~ 10"^ 


0.65 ±0.2 


2.25 ±0.55 


~ 10"* 0.15 


1.0 


0.95/19 


S5 0716+714 


30-31 October, 2000" 


1.25 ±0.08 


16.4 ± 1.2 






1.1 


2.9 


1.99/76 


30-31 October, 2000" 


2.40i-i 


10.2 ±0.45 


^•"'-'-0.35 


"•^-3.0 


0.7 0.4 


3.3 


0.92/77 



Table 2. X-ray spectral parameters and fluxes. *: cts keV ^ cm ^ s ^. **: erg cm ^ s ^. The values are obtained performing 
the fit with the absorption parameter fixed at the galactic value N^i — 1.18 x 10^" cm~-^ for OQ 530 and A^h — 3.9 x 20^" cm^^ 
for S5 0716+714 (Dickey & Lockman, 1990). fit performed only on the LECS+MECS data, ^ added also flie PDS data to the 
fit. The reported uncertainties are the 90% confidence ranges for one parameter of interest, i.e. Ax^ = 2.7. 




0.2 0.5 1 2 5 10 1 10 

channel energy (keV) channel energy (keV) 



Fig. 1. Left panel: LECS+MECS S5 0716+714 spectrum fitted by a single power law model. Positive residuals are evident below 
^ 1 keV; they get larger towards lower energies. Right panel: LECS+MECS+PDS S5 0716+714 spectrum fitted by a two power 
law model. In both cases we fitted the data with the absorption parameter fixed at the galactic value TVh = 3.9 x 10^" cm^^ 
(Dickey & Lockman, 1990). 



tra with single power laws, keeping the absorption parameter 
fixed at the galactic value: N^i = 3.9 x IQ-^" cm^^ (Dickey 
& Lockman, 1990). The LECS+MECS spectrum is badly fit- 
ted by this model {xi/d.o.f. = 1.99/76), with large positive 
residuals below 0.5 keV (see Fig. Therefore we fitted it 
again with a two power law model. This model fits the data 
much better (ai = 2.5 aa = 0.7, xl/d.o.f. = 0.86/74). 
The two power laws cross at ~ 1 keV. The MECS+PDS spec- 



trum is instead fitted quite well by a single power law model, 
with a spectral index somewhat steeper than the second spec- 
tral index of the LECS+MECS two power law fit {a 0.9, 
Xr/d-O-f- — 1.08/36). However, we note that due to the weak 
signal in the PDS, the fit is clearly dominated by the MECS 
data. A detailed analysis of the PDS light curve in search for 
spikes gave negative results. This, together with the lack of 
other obvious X-ray sources inside our field of view, confirm 
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0716 + 714 




Log i/(Hz) 

Fig. 2. X-ray spectral energy distribution of S5 0716+714 dur- 
ing the October 2000 BeppoSAX observation (filled circles). 
We plotted also November 1996 (open triangles), November 
1998 (open squares) BeppoSAX (Giommi et al. 1999) and 1991 
ROSAT SEDs (Cappi et al. 1994). 



variation could be even larger since our measurements stopped 
before detecting any fading. The source displayed a similar be- 
haviour also in the 1.5-3.0 keV MECS band (not plotted here). 
In spite of the formal break that seems to be at around 1 .5 keV, 
the emission in the two softer energy bands are clearly con- 
nected. Clearly, the soft X-ray component component is still 
relevant in the 1-3 keV band. However, by selecting only the 
counts above 2 keV we find that the variability is much less 
pronounced. The variability is essentially absent in the 3-10 
keV band, where the hard X-ray component is dominant (see 
Fig.|3}- A similar behavior was already detected in a previous 
BeppoSAX observation (Giommi et al. 1999). 

3. Optical, Radio and TeV data 

3.1. OQ530 

The optical Rc light curve since February 1996 of the source 
is shown in Fig.|4] It is characterized by a large variability 
with fast and strong flares (of the order of half magnitude in 
a few days). The average trend seems to be decreasing un- 
til the beginning of 1999, when an increasing trend seems to 
begin. The average value was Rc =15.3 in the period Jan 
1997 - Aug 1999 (JD 400 - 1400) so that our observation 
at Rc =14.3 on Feb 14, 2000 was clearly indicative of a re- 
markably bright state. During the first BeppoSAX pointing on 
March 3 the source had already dropped at i?c'= 14.9 and re- 



us on the reality of the PDS detection. We then fitted the whole 
LECS+MECS+PDS spectrum with a two power law model, 
finding a good fit to flie data (xl/d-o.f. = 0.92/77). The two 
power laws cross at 1.5 keV. In Table 2 we report the best-fit 
parameters. 

Although S5 0716+714 was not at its maximum recorded 
optical luminosity (Rc ~12.5), measured a few days before 
our -X-ray observation, it was still about 0.8 mag brighter than 
during the two previous BeppoSAX exposures (Giommi et al. 
1999). It was brighter than the previous observations also in 
the X-ray band (see Table|3l During our observation the source 
was in an X-ray state higher in flux but similar in shape to that 
detected by ROSAT in 1991, when it was revealing again a con- 
cave spectrum with an energy break at ~ 0.8 keV (Cappi et al. 
1994, Comastri et al. 1995). The historical SED in the X-ray 
band is displayed in Fig. |2] In all observations the transition 
between the steep soft X-ray component and the much harder 
spectrum at higher X-ray energies is clearly detected. During 
our observation the soft X-ray spectrum is steeper (ai = 2.4), 
than during the previous BeppoSAX observations (a — 1.7 in 
1996; a = 1.3 in 1998). This seems to confirm the trend no- 
ticed by Giommi et al. (1999) of a steeper soft X-ray spectrum 
when the source is brighter. 

We performed the temporal analysis of our observation in 
file three different energy bands: 0.15-1.0, 1.0-3.0 and 3.0-10 
keV. In Fig.|3]we display the two LECS and the higher en- 
ergy MECS light curves. They seem almost constant for the 
first 20 hours of the observation. Then the two LECS fluxes 
increase by a factor of three on a time scale of one hour. This 



S5 0716+714 (30-31/10/2000) (Bin time=600s) 
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Fig. 3. LECS 0. 15-1 keV (top panel), 1-3 keV (mid panel) and 
MECS 3-10 keV (bottom panel) light curves of S5 0716+714 
during the BeppoSAX observation of October 2000. The time 
bin is of 10 minutes. Bins with less than 20% effective expo- 
sure time have been excluded. Time axis is hours starting from 
October 30, 2002. 
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Fig. 4. The optical, Rc magnitude, light curve of OQ 530 since 
February 1996. Open triangles are data from the Perugia obser- 
vations (the typical error for these measurements is of the or- 
der of 0. 1 magnitude), filled circles are from Roma (the typical 
error for these measurements is of the order of 0.01-0.02 mag- 
nitude). The two arrows indicate the dates of the BeppoSAX 
pointings, activated by a flare-like episode, but carried out 
when the source had already faded to Rc ~ 15 mag (see text). 

mained stable within 0.03 mag; during the second one (March 
26) the source was about at the same optical luminosity but 
we could observe it only for a short time due to bad weather 
conditions. Between the two BeppoSAX pointings the source 
remained substantially constant, with small (^0. 1 mag) varia- 
tions. The energy spectral slope (F^ oc i/~") from our BVRI 
observations was 1.63±0.06, a typical value for the source at 
this luminosity level (D' Amicis et al. 2002). On the other hand, 
the spectral shape (and flux levels) is affected by the contribu- 
tion of the host galaxy, and in Fig. 6 the simultaneous optical 
points have been galaxy-subtracted, using the host magnitude 
given in the Introduction. 

Very Long Baseline Interferometric radio observations with 
EVN at 5 GHz were performed one year before (February 

1999) and one year after (June 5, 2001) the BeppoSAX point- 
ing by some of us (Massaro et al. 2002). The source showed 
a core-jet structure, already detected in 1996 (Fomalont et al. 

2000) : the core had a flux density of 440 mJy at that epoch, 
but dropped to 204 mJy on February 1999 and rose again up to 
685 mJy on June 2001. This increase, larger than a factor of 3, 
of the core radio brightness supports the development of strong 
activity following the optical burst of February 2000. 

OQ530 has been pointed in the TeV band by the HEGRA 
telescopes from July 1 to July 16 1999, and on March 5 2000. 
The total exposure time is of 9.7 hours, a quarter of which is 
from the March observation (1 day after the BeppoSAX obser- 
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TIME 

Fig.5. The 1.5 years light curve of S5 0716+714 at three radio 
frequencies: triangles are data at 4.8 GHz, open circles data at 
8.0 GHz, and x denotes data at 14.5 GHz. The measurements 
are given as daily averages. The source show erratic variability 
during all year 2000, then there seem to be a steady decline in 
the fluxes. 

vation). The source was never detected with a total upper limit 
UL(99%, E> 1.1 TeV) = 1.35 x 10"i2 photons cm^^ g-i 
(Aharonian et al., HEGRA collaboration, 2002, in preparation). 

3.2. S50716+714 

As stated in the Introduction, the source was very bright the 
week before the BeppoSAX observation: Rc — 12.55 ± 0.02 
mag on October 22 and 12.58±0.03 on October 28. Only few 
optical photometric data, taken during the X-ray pointing, are 
available, but these are enough to derive the SED and to com- 
pare it with that of previous observations. In particular, Rc 
and Ic were measured equal to 13.05±0.04 and 12.55±0.04 
mag at 22:40 UT (October 30) with the 0.5m reflector of the 
Astronomical Station of Vallinfreda (Rome). On November 1 
22:10 UT its luminosity has faded at i?c=13.31 (VaUinfreda), 
indicating that over the entire period S5 0716+714 showed 
brightness variations on intraday time scale with amplitude of 
a few tenths of magnitude, as is usually observed in this source 
(Nesci et al. 2002). 

During the BeppoSAX observation the radio flux of 
S5 0716+714 was monitored at the UMRAO, using the proce- 
dures and calibrators described in Aller et al. (1985). At 14.5 
GHz, from October 28 to November 1, the flux density re- 
mained constant with a mean value of 1 . 1 Jy. Other measure- 
ments at 8.0 GHz, performed from October 1 to 21 give a flux 
density in the range 0.60 — 0.85 ± 0.1 Jy. In Fig.|5lwe plot the 
1 .5 years radio light curve of S5 07 16+7 14 at three frequencies. 
The source seems to be in an active state up to the end of year 
2000, then there seems to be a steady decline in the fluxes. 

During our BeppoSAX observation S5 0716+714 has been 
observed also in the TeV band by the HEGRA telescopes for 
about 2 hours on October 31 2000. However, the source has 
not been detected with an upper limit UL(99%, E > 1.7 TeV) 
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OQ530 



Date 


Mission 




a2 




19 _9 _1 

xlO erg cm ^ s ^ 


Ref. 


December 1980 


Einstein 


o.48li:!; 




0.21 




1 


July 1990 


ROSAT 


1.04 ± 0.05 




0.32 




2 


February 1999 


BeppoSAX 


55+"-^'' 




0.09 


0.67 


3 


3-4 March 2000 


BeppoSAX 


0.55 ± 0.2 




0.14 


1.1 





26-27 March 2000 BeppoSAX 0.75 ± 0.2 (U6 03 



S5 0716+714 



8 March 1991 


ROSAT 


2.67 ±0.12 


1.00 


0.78 




4 


14 November 1996 


BeppoSAX 


1.7 ±0.3 


0.96 ±0.15 


0.13 


1.4 


5 


7 November 1998 


BeppoSAX 


1.3 ±0.4 


0.73 ±0.18 


0.19 


2.6 


5 


30-31 October 2000 


BeppoSAX 


2A0t\\i 




1.1 


3.3 





Table 3. OQ530 and S5 0716+714 X-ray spectral parameters. (1) Worrall & Wilkes (1990); (2) Comastri et al., 1995; (3) 
Padovani et al., in prep.; (4) Cappi et al. (1994); (5) Giommi et al. (1999). 



= 2.5 X 10~^^ photons cm~^ (Aharonian et al., HEGRA 
coUaboration, 2002, in preparation). 

4. The Spectral Energy Distributions 

We have used a homogeneous, one-zone synchrotron inverse 
Compton model to reproduce the SEDs of our sources. The 
model is very similar to the one described in detail in Spada et 
al. (2001), it is the "one-zone" version of it. Further details can 
be found in Ghiselhni, Celotti & Costamante (2002), where the 
same model has been applied to less powerful BL Lacs. The 
main assumptions of the model are: 

- The source is cylindrical, of radius R and thickness Ai?' = 
R/T (in the comoving frame, where T is the bulk Lorentz 
factor). 

- The source is assumed to emit an intrinsic luminosity L' 
and to be observed at a viewing angle 6 with respect to the 
jet axis. 

- The particle distribution N{j) is assumed to have the slope 
n [N{'y) oc 7~"] above the random Lorentz factor jc, 
for which the radiative (synchrotron and inverse Compton) 
cooling time equals AR'/c. The motivation behind this 
choice is the assumption that relativistic particles are in- 
jected in the emitting volume for a finite time, which we 
take roughly equal to the light crossing time of the shell. 
This crossing time is roughly equal to the time needed for 
two shells to cross, if they have bulk Lorentz factor dif- 
fering by a factor around two. The electron distribution is 
assumed to cut-off abruptly at 7max > 7c- We then as- 
sume that between some 7min and 7c the particle distri- 
bution AT (7) oc 7~("~i). This choice corresponds to the 
case in which the injected particle distribution is a power 
law (oc 7~("~i)) between With7min < 7c- 
Below 7niin we assume N{j) oc 7^^. The value of 7max 
is not crucial, due to the fact that the electron distribution, 
for our sources, is steep. It has been chosen to be a factor 
^ 100 larger than 7niin- 



The simultaneous optical observations help in defining the 
low frequency peak of the SED, which in both sources must 
be located at frequencies lower than the optical. This helps to 
constrain the values of the input parameters of our model. Since 
we determine 7pcak as the energy of those electrons that can 
cool in the injection time (i.e. 7peak = 7c)> there is a relation 
between fpeak and the cooUng rate. If the latter is dominated by 
synchrotron emission, the location of t'pcak constrains the value 
of the magnetic field. This in turn fixes the amount of inverse 
Compton radiation. According to our simultaneous data, the 
slope of the optical emission for both sources is steep (i.e. a > 
1), requiring the synchrotron part of the spectrum to peak below 
the optical band. On the other hand, the infrared data, even if 
not simultaneous, suggest that i^poak cannot be far below the 
optical. This motivates the choice of the adopted i^peak as listed 
in Table 4. 

The remaining degree of freedom for the choice of the in- 
put parameters is due to the value of the beaming factor (i.e. 
r and the viewing angle 9), and to the redshift. Another im- 
portant input parameter is the size of the emitting region. The 
extremely rapid variability shown by S5 0716+715 during our 
observations refers to the soft X-ray flux, while timescales of 
the order of a ~day are typical in the optical (see e.g. Ghisellini 
et al. 1999), at frequencies closer to the synchrotron peak. The 
size of the emitting region is therefore constrained to be less 
than one light-hour or a light-day (i.e. R < ct^a.rS/ (1 + z)). 
Although a one-zone homogeneous model is forced to use the 
minimum variabihty timescale observed at any band to con- 
strain the size, it is also clear that this model is a simplifica- 
tion of a scenario which may be more complex. Our choice of 
R = 2 X 10^^ and 5 ~ 17 corresponds to a minimum vari- 
ability timescale of tyar ~ 11(1 + 2;) hours. Assuming a more 
compact homogeneous source, while more in agreement with 
the soft-ray variability, implies a more dominant self Compton 
component, unless a larger beaming factor is also assumed. For 
the redshift, since only the lower limit z = 0.3 is known for 
S5 0716+714, we have assumed z = 0.3 in all our calculations 
for S5 0716+714 but one case, for which z = 1 has been as- 
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sumed, in order to check what are the output parameters mostly 
affected by the choice of a particular redshift. 

We have tried to reproduce the three SEDs of S5 0716+714 
by changing the minimum possible number of input parame- 
ters: the most significant change is in the injected luminosity, 
which in Oct 2000 is found to be twice as much as in 1996 and 
1998. 

An important consequence for our modeling is that the SSC 
component alone cannot reproduce the observed emission in 
the EGRET band in S5 0716+714, as shown in the bottom 
panel of Fig. 7. Since these 7-ray data are not simultaneous, 
there is in principle no need to reproduce them. However, since 
the source should be in a bright state, it may be likely that the 
level of 7-ray emission during our campaign was similar to 
what EGRET previously found (Lin et al. 1995). If we want 
to account for a flux in the 0.1-10 GeV band similar to what 
found by EGRET, we should then consider the possibility that 
the high energy emission can be produced by the scattering of 
external soft photons, either produced in a putative broad line 
region or in the accretion flow. We have allowed for some ex- 
ternal Compton radiation by assuming the presence of external 
photons produced by, e.g. a broad line region of relatively low 
luminosity, as listed in Table 4. The effect of this mechanism is 
to produce the second "hump" at the largest 7-ray energies, as 
shown by the solid line in Fig. 6 (for OQ 530) and by the top 
panel of Fig. 0(for S5 0716+714). This high energy bump is 
due to the fact that the seed photons due to the emission lines, in 
the comoving frame, are at the typical frequency of lO^^F Hz, 
larger than the synchrotron peak frequency in the same frame, 

which is t^pcak/<^ ^ t'pcak/r). 

The possibility of a broad line region in objects such as 
S5 0716+714 is questionable: the line emission could be vari- 
able (as in BL Lac itself), or always very weak (as seems to 
be the case for low power, high frequency peaked BL Lacs). 
On the other hand, to produce extra inverse Compton emission, 
other processes could be important, as in the "mirror model" 
discussed by Ghisellini & Madau (1996). Therefore the values 
listed in Table|4] which formally refer to the broad Une region, 
could instead be appropriate for other possible contributors of 
seed photons. Note that for OQ530 there are only upper limits 
for the 7-ray flux, which are not stringent enough to constrain 
the presence of any external component (compare the solid and 
dashed lines). 

The input parameters are listed in Table |4] and the models 
are shown in Fig.|6land Fig.0 The applied model is aimed at 
reproducing the spectrum originating in a limited part of the 
jet, thought to be responsible for most of the emission. This 
region is necessarily compact, since it must account for the fast 
variability shown by all blazars, especially at high frequencies. 
The radio emission from this compact region is strongly self- 
absorbed, and thus the model cannot account for the observed 
radio flux. This explains why the radio data are systematically 
above the model fits in the figures. 

5. Discussion 

We presented the results of BeppoSAX ToO observations of 
two BL Lac objects observed while they were in a high state. 




10 15 20 25 

Log u [Hz] 



Fig. 6. The SED of OQ530. Star symbols are from the NASA 
Extragalagtic Database (NED). The solid line refers to the SSC 
plus EC model, while the dashed line corresponds to a pure 
SSC models, with no contributions from seed photons gener- 
ated externally to the jet. All other parameters of the two mod- 
els are the same. 



Once again these observations were triggered from optical 
monitoring. In our ToO program we have probably been biased 
toward sources that have an higher optical variability. These 
should be the blazars that have the synchrotron peak in the IR- 
optical band and these are of course essentially LBL or inter- 
mediate blazars. In particular the two sources studied in this 
paper belong to the intermediate BL Lac class in terms of their 
spectral energy distribution. In this regard they resembles ON 
231 (Tagliaferri et al. 2000): all these three BL Lacs show ev- 
idence for a concave X-ray spectrum in the 0.1-10 keV band, 
a signature of the presence of both the steep tail of the syn- 
chrotron emission and the flat part of the inverse Compton 
spectrum. Another source that we studied more than once and 
in which we detected both components, with a concave X- 
ray spectrum, or only the Compton component, is BL Lac it- 
self (Ravasio et al. 2002). A common result that we found in 
ON 231, BLLac and S5 0716+714 is that fast variability was 
detected for all the three sources only in the synchrotron com- 
ponent. This can be interpreted with the presence in the X-ray 
band of a Compton component (slowly variable on time scale 
of months), and the tail of a synchrotron component with fast 
and the erratic variability. 

The Compton emission we see in the X-ray band is well 
below the Compton peak and it is produced by low energy 
electrons scattering low frequency synchrotron photons. The 
short timescale variability seen in the synchrotron part can be 
reproduced by changing the slope of the injected electron dis- 
tribution, without affecting the total injected power. This is in 
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Table 4. Model Input Parameters. Column (1): name of the source; column (2): observation date; column (3): intrinsic luminosity 
L'\ column (4): luminosity of broad emission lines Lblr: the first value (i.e. zero) corresponds to assume no contribution 
of external photons to the formation of the IC spectrum; column (5): dimension of the broad line region i?BLR; column (6): 
magnetic field B; column (7): size of the emitting region R\ column (8): bulk Lorentz factor F; column (9): viewing angle (in 
degrees); column (10): beaming factor 6; column (11): slope of the particle distribution n; column (12): minimum Lorentz factor 
of the injected electrons 7inin; column (13): Lorentz factor of the electron emitting at the peaks, 7pcak; column (14): synchrotron 
peak frequency t'pcak- Note that 7pcak and v-pcak are derived quantities and not input parameters. 



Name 


date 


L' 


isLR 




B 


R 


r 


e 


5 


n 


7min 


Tpcak 


^pcak 






ergs"^ 


erg 


cm 


G 


cm 














Hz 


0716 (z=0.3) 


Oct 2000 


5.0e42 


0-5e42 


3el7 


2.5 


2el6 


15 


3.4 


16.7 


3.7 


5.0e2 


1.8e3 


5.1el4 


0716 (z=0.3) 


Nov 1998 


2.7e42 


0-5e42 


3el7 


2.5 


2el6 


15 


3.4 


16.7 


3.9 


4.0e2 


2.0e3 


6.2el4 


0716 (2=0.3) 


Nov 1996 


2.2e42 


0-5e42 


3el7 


3.0 


2el6 


15 


3.4 


16.7 


3.9 


5.0e2 


1.5e3 


4.2el4 


0716 (z=l) 


Oct 2000 


4.0e43 






5 


2el6 


15 


3 


18.5 


3.7 


1.2e3 


1.2e3 


5.0el4 


OQ530 


Mar 2000 


8.0e41 


0-5e42 


2el7 


2.8 


lel6 


10 


5.0 


11.4 


3.9 


3.3e2 


2.2e3 


5.8el4 



line also with the fact that we do not see large shifts of the 
synchrotron peak frequency during our observations, that are 
usually performed when the source is in a higher state than the 
previous observations (e.g. Fig.^l see Tab. 4). 

In the framework of the blazar sequence proposed by 
Fossati et al. (1998) the two sources studied here have bolomet- 
ric luminosities intermediate between HBLs (with bolometric 
luminosities L ~ 10"*^ erg s^^) and powerful broad line blazars 
(with bolometric luminosities up to L ^ 10^^ — 10^^ erg s^^, 
see Fig. 12 in Fossati et al. 1998). In the case of S5 0716+714, 
assuming z — 0.3 and a cosmology with Hq — 65 km s^^ 
Mpc~^, = 0.3 and J7a = 0.7, the bolometric luminos- 
ity is - (3 - 5) X 10^*^ erg s~\ while in the case of OQ 530 
we have (2 — 10) x 10"*^ erg s~^. The values of their 7poak 
and the comoving radiation plus magnetic energy density lo- 
cate our objects along the correlation found by Ghisellini et al. 
(1998) and GhiselUni, Celotti & Costamante (2002), as shown 
in Fig. |8l Notice also that the two objects are located at the 
joining point between the two branches (i.e 7pcak « f7~^ and 
7pcak oc U^^/'^) proposed by Ghisellini, Celotti & Costamante 
(2002) as characterizing all blazars (we also report for com- 
parison one case with z = 1 for S5 0716+714). Note also that 
the 3 states of S5 0716+714 studied here correspond to this ob- 
ject "moving" along the 7pcak-C^ correlation, even if the small 
dynamic range of the variations makes this result tentative at 
best. 

Our results show the importance of monitoring these vari- 
able sources and of performing multiwavelength observations, 
including a X-ray energy band as wide as possible, while they 
are either in a high or in weak state, to obtain the largest in- 
formation on the dynamical evolution of the emission compo- 
nents. 
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applied to the three labeled observing campaigns. For all mod- 
els z = 0.3 was assumed. The bottom panel shows pure SSC 
models, with no contributions from seed photons generated ex- 
ternally to the jet. The other input parameters are the same 
(solid lines), while z = 1 is assumed for the dashed line. The 
square symbols refer to not simultaneous observations taken 
from NED and from Lin et al. (1995). 
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is the sum of the radiation and magnetic energy densities, as 
measured in the comoving frame. Note that our two sources 
are located in the region where the two branches (proposed by 
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report also one case with z — 1). Adapted from Ghisellini, 
Celotti & Costamante (2002). 
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